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Abstract 

A many-body theory based on nonequilibrium Green functions, in which transport and optics are 
treated on a microscopic quantum mechanical basis, is used to compute gain and absorption in the 
optical and THz regimes in quantum cascade laser structures. The relative importance of Coulomb 
interactions for different intersubband transitions depends strongly on the spatial overlap of the 
wavefunctions and the specific nonequilibrium populations within the subbands. The magnitude 
of the Coulomb effects can be controlled by changing the operation bias. 
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I. INTRODUCTION 



The quantum cascade laser (QCL) has evolved, since its first experimental realization 
[lj], into an important device for infrared applications. Recently, lasing in the THz regime 
has been demonstrated, thus opening further perspectives P . The operation of these semi- 
conductor heterostructure devices is based on optical intersubband transitions within the 
conduction band, which are the focus of this work. The importance of many-body effects 
on interband transitions is a well established fact |3|. The intersubband case has also been 
studied both experimentally and theoretically [4] and is a topic of current research P, |(| . 
Typically, idealised cases of isolated quantum wells with equilibrium distributions for the 
carriers are considered in a two subband scenario, to facilitate the computational challenges. 
Assuming that only one level is occupied, exchange and depolarisation effects were found 
to be very relevant for both gain and absorption [2J. Recent work has demonstrated the 
importance of fullv solving the integral equation for the optical susceptibility for quantum 
wells 0, 0, 0, 0, 10]. While standard semiconductor lasers based on interband transitions 
can be modelled by equilibrium distributions for electrons and holes with different quasi- 
Fermi levels h|, this is not appropriate for QCLs where a full kinetic theory is necessary 
to correctly determine the nonequilibrium populations of the multiple level system. 

The main issues of our paper are: (i) the question how these many-body effects are 
affected by the inherent nonequilibrium distribution in a QCL, where the interplay of a 
large number of different levels is crucial for the operation of the device; (ii) to what extent 
many-particle corrections can be modified in a given structure out of equilibrium and in a 
complex wavefunction scenario. 

We show that a characteristic QCL design provides rather small overall Coulomb cor- 
rections to the gain transitions in the infrared, while significant level shifts occur in the 
THz absorption region. The latter depend strongly on the bias applied to the structure, 
which opens a possibility for a systematic variation of Coulomb effects in an experimental 
study of a complex nonequilibrium system. Here a wider variation of parameters is possible 
compared with previous experimental studies in quantum wells by either applying a bias 



comparea witn previous experimental sti 
[3, [3] or optically pumping a laser [l^ . 



Our findings thus give evidence of the need of our more advanced approach to model 
complex nonequilibrium studies. The study of simple structures in ref. 0] would suggest 
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similar strengths of Coulomb corrections in both gain and absorption regimes. Although 
those results are qualitatively correct for the cases previously studied, they can not be used 
for more complex structures of current interest, as demonstrated in our analysis, since the 
Coulomb matrix elements and occupation functions are extremely different for gain and 
absorption transitions. That large discrepancy is explained here for the first time to the 
best of our knowledge. 

II. MAIN EQUATIONS 

Our calculations are based on the stationary state under operating conditions which is 
obtained by a self-consistent solution of the quantum kinetic equations Q| • Here scattering 
with phonons and impurities are included via the respective self-energies in the Dyson equa- 
tions, and electron-electron interaction is treated in mean field approximation. The lattice 
temperature only enters via the phonon bath at 77 K, which is assumed to be in thermal 
equilibrium, neglecting hot phonon effects [18]. For each subband, the nonequilibrium occu- 
pation functions, f u (k) are calculated from the diagonal terms in the carrier Green functions 
G^(k, uj) which are expanded in terms of eigenstates <fi v {z) of the single particle Hamilto- 
nian including the mean field. (These states will be referred to as levels in the following.) 
The calculated dephasing 2Y U is the FWHM of the spectral function, ^s{G r ^(k = 0, u)} and 
E u (0) the corresponding location of the peak of the spectral function so that the electronic 
states are well defined quasi-particles, whose properties are evaluated on a microscopic basis. 
The optical absorption a at a given photon energy hut can be calculated from the imaginary 
part of the optical susceptibility, 



Here n b denotes the background refractive index, c is the speed of light, p vlx is the transition 
dipole moment between the subbands v and fi, which are labelled v — 1, 2, . . . with increas- 
ing energy. In our calculation we select the 9 lowest subbands per period. The definition of 
a period (and thus the choice of level 1) is of course arbitrary in such structures and we must 
also take into account transitions between neighbouring periods. Numbering the levels in 
the next period to the left by 10-18, we consider all transitions (1, . . . 9) «-> (1, . . . 18) in the 
computations. We evaluate the susceptibility function, Xu^(k,u), by the carriers' Keldysh 




(1) 
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Green's function, G, whose time evolution is described by a Dyson equation, with Coulomb 
interactions as well as other scattering mechanisms included in self-energies, E The 
general theory includes carrier-carrier scattering, diagonal and non-diagonal dephasing. In 
the numerical results presented here, we keep only the subband shift, depolarisation and 
exchange terms in the electron-electron interaction (Hartree-Fock). Similarly to the inter- 
band case, but with additional terms due to depolarisation as well as full consideration 
of a Hartree (mean field) term in the quasi-free particle term ^(k), the resulting integral 
equation for the susceptibility function reads, 



- (f v (k) - / M (k)) Y,x^',")v( 



v v fx n 

k-k' 





k' 

p„ M (/„(k) - /„(k)) . (2) 



kVk 

Terms with more then two different indices vanish identically in the limit of isolated wells 
with infinite confinement barriers and remain small in the more complex structure studied 
here. The bare Coulomb interaction is given by 

ffiu\(3\ f 27ce 2 exp(-\k-k'\\z-z'\) 

V k - k/ / = J M^M 2 ) e A\k-k'\ ( 3 ) 

where we have introduced the normalisation area A, and the background dielectric constant 
cq. The explicit expressions actually programmed here are given in the Appendix. 



e^(k) =E u (k) - E,(k) - AOO V ( V ™£) + E ^( k ') V ( 



k-^ 



k-k' 



(4) 



k' 

is the energy difference between the levels renormalised by the exchange interaction, which 
we refer to as subband shift in the following. The second term on the left-hand side of 
Eq. (J2J) gives rise to the depolarisation shift, while the last term (exchange contribution) is 
analogous to the excitonic coupling term in interband transitions The equations above 
reduce to those of Ref. [19| in the two-band limit at equilibrium. Note that by running over 
all subband indices v ^ \i in Eqs. (jll2|) . we do not apply the rotating wave approximation. 
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III. NUMERICAL RESULTS AND DISCUSSION 



For numerical applications of the theory, we choose the very successful prototype intro- 
duced by Sirtori et al. see Fig. [TJ From a systematic study of doping dependence 1(| 
it was found that the performance is best for a doping around 6.4 x 10 11 /cm 2 which is used 
throughout this work together with a lattice temperature of 77 K unless stated otherwise. 

In Fig. El we compare the absorption spectra evaluated with and without Coulomb effects 
for different biases. We find that gain sets in around 0.2 V/period at energies %u ~ 130 
meV (infrared region) in good agreement with the experimental findings 16j. In addition 
there is a pronounced absorption around 30 meV (THz- region). While the Coulomb effects 
are hardly visible in the gain region, significant shifts of peak height and position occur in 
the THz absorption lines, in particular at low bias. 

In the following we present a systematic study of the Coulomb terms in Eq. (J2J), which 
crucially depend on the Coulomb matrix elements and the occupation differences, which 
are displayed in Table H] for the dominant transitions. A glance at the values n v — 
[which roughly scale with f u (k) — / M (k)] shows that the structure is out of equilibrium, as 
necessary for laser operation. In order to quantify changes in Coulomb effects, we express 
the matrix element for the depolarisation term as Lf, ep = 9 ^° 2 A , so that it reduces to 
the well width in the limiting case of a quantum well with infinite confinement barriers for 
transitions between the two lowest subbands. This reflects the common knowledge that the 
depolarisation shift increases with well width, i.e. increasing spatial extension of the wave 
functions. For a QCL the scenario is much more involved, as the spatial overlap of the wave 
functions essentially contributes to the Coulomb effects as well. Similar quantities can also 
be defined for the exchange and subband shift, namely LfZ, and L\f u . Explicit expressions 
are derived in the Appendix and numerical values are given in Table |l] In contrast to 
L dep , both L ex , and L b9 are inversely proportional to the corresponding Coulomb matrix 
elements. We can then expect that relatively large L dep and small L ex should correspond 
to large dipole moments, or in other words, to strong overlap between the wavefunctions 
of the levels involved in the transition. Only the wavefunction for level v appears in L\f v . 
Thus, its value remains within the same order of magnitude irrespective of the transition. 
Note that the significant difference between the values of L dep , L ex , and L bg as well as for the 
different transitions further demonstrates that results from a simplified two-level system in 
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a quantum well, where these lengths are identical by definition, cannot be easily generalised 
to more complex scenarios, like the real QCL structure discussed here. 

In the THz region, and the dipole moment p UfJi tend to decrease with increasing bias, 
due to the changes in the wave function overlap. Furthermore the population decreases in 
the injector. This explains the general trend that the Coulomb effects become weaker with 
increasing bias. At 0.2 V/period, the concentration of the population in the lowest injector 
level 7 leads to a very strong absorption feature in the THz region although L^f is smaller 
than at 0.12 V/period. In the infrared region, the larger bias increases the occupation of the 
subbands which are responsible for laser action. Furthermore both Lf,^ and p Ufl increase 
with bias as the injector levels leak into the active region. Nevertheless L d ^ is much smaller 
than in the THz absorption region, and furthermore the compensation between the different 
many-particle effects (see below) in the dominating gain transitions leads to very small 
resulting shifts in Fig. El The contrast between the interesting Coulomb shifts in the low 
energy side and the almost vanishing shift in the high energy gain side is dramatic. This is the 
reason why gain calculations without Coulomb corrections [12f could explain experimental 
gain spectra in this structure reasonably well, but will fail in the absorption region. The 
main contributions to the absorption in the THz region are in fact transitions involving the 
two lowest injector levels at each bias. In other words, levels 6 and 7 for 0.12 V/period, 
levels 7 and 8 for 0.2 V/period and levels 7 and 9 for 0.3 V/period. 

The interplay between the different Coulomb corrections in the THz region is analysed 
in more detail in Fig. El We select the dominant transitions for each bias to highlight the 
different contributions more clearly. The general tendency is that the exchange contribution 
imparts a red-shift while subband shift and depolarisation terms lead to blue-shifts. The 
depolarisation shift decreases with bias (in accordance with the discussion above), while the 
tendency of the exchange contribution and the subband shift is less clear. These two terms 
almost cancel each other, so that the blue-shift due to depolarisation dominates the be- 
haviour. Because the subband shift slightly overcompensates the exchange contribution, the 
full Hartree-Fock curve exhibits a small additional blue-shift, in particular for 0.2 V/period. 

Figure E] shows details for the two dominating gain transitions at 0.3 V/period. We no- 
tice first that all many-particle effects reverse sign for gain transitions. This is due to the 
sign change in the occupation differences entering Eqs. (I2I4D |7|. Both transitions exhibit a 
red-shift due to all the Hartree-Fock corrections. The shift is stronger for the (6,7) transition 
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although the depolarisation length is much smaller compared to the (6,8) transition, see Ta- 
ble HI For gain transitions a significant part of the upper laser level is located in the injector, 
resulting in a weaker spatial overlap of the wave functions M and <f> v . Thus the Coulomb 
matrix elements entering exchange and depolarisation corrections (such as L^) are reduced 
compared to the THz absorption within the injector region. In contrast Coulomb matrix 
elements of the type V(fi, /i, //, //) enter the subband shift (JH), which therefore dominates 
the behaviour in the gain region. Indeed, a lookup at Table |l] shows that L dep is small, L ex 
is large, consistently with the small dipole transition moment, while L bg is relatively small, 
leading to the scenario depicted in Fig. HJ 

The spectra of typical QCLs are found to result from a delicate interplay of different 
transitions with strongly varying transition parameters. Standard approaches for idealised 
quantum wells cannot describe this situation, where the full nonequilibrium distribution ob- 
tained from quantum transport or possibly semi-classical Monte-Carlo [20] simulations, 
must be combined with the correct wavefunctions for a consistent description and under- 
standing of the detailed spectra. Table |U shows that L^ff and p Ufl can be easily modified 
within the same structure by changing the bias. This allows for a systematic experimental 
study of the magnitude of Coulomb effects as a function of transition parameters. Thus QCLs 
constitute promising structures to investigate the influence of Coulomb effects on intersub- 
band transitions. A central issue of such experimental studies as well as future theoretical 
work will be the importance of higher order Coulomb corrections together with a more strin- 
gent treatment of screening compared to our Hartree-Fock approach. Recent work indicates 
that compensations between diagonal and nondiagonal terms may strongly reduce the 
importance of higher-order electron-electron scattering terms in simple model structures. 
This indicates that neglecting scattering in our Hartree-Fock approach may be a reasonable 
approximation, although this requires further detailed investigation. 

Furthermore we illustrate the importance of many body effects as compared to the free 
carrier model as a function of both doping densities and temperatures. In order to keep our 
discussion as close as possible to experimental realizations, we chose densities and temper- 
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atures compatible with the samples described in Refs. 

Figure |S] shows the absorption spectra with and without Hartree-Fock Coulomb correc- 
tions for different doping densities. Our free carrier model includes the mean field of the 
ionised dopants and as the density increases, they shift the peak of the main absorption 
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transition to the blue. The many-body effects impart a further intersubband blue shift and 
increase in oscillator strength of the main transitions. Although both the subband shift and 
overall oscillator strength increase with carrier density, the impact of the Coulomb correc- 
tions is similar for all densities and follows our detailed analysis of the N=6.4 x 10 11 /cm 2 
case given above. 

The density dependence study is complemented with Fig. |B] that depicts the evolution 
of gain (-absorption) spectra with density. Maximum gain is reached around N=4.6 x 
10 n /cm 2 . The lower panels show that Coulomb corrections are larger for the larger density 
and the inset depicts the evolution of peak gain position with doping. Although the cor- 
rections are small it is interesting to note that Coulomb corrections lead to smaller increase 
in peak gain position on the high density, making the results slightly closer, qualitatively, 
to the experimental findings of Ref. [lfj]. However, Coulomb effects only are not enough to 
explain their results and a more detailed analysis is required. 

We close our analysis with Fig. [7| where we show the influence of temperature on the 
calculated curves with (solid) and without (dot-dashed) many-body effects. The top and 
lower panels are respectively for phonon bath temperatures T=77 and 180 K. In both panels, 
the bias is V=0.2 V/period and the carrier density N= 6.4 x 10 11 /cm 2 . 

The rise in temperature increases the spectral broadening and promotes a redistribution of 
population difference, favouring the gain transitions As a consequence, the resulting subband 
shift is not so pronounced in the low energy side in comparison to T=77K. However, since 
different transitions have different Coulomb corrections, the inclusion of Coulomb effects at 
T=180 K affects the qualitative shape of the spectrum more strongly at higher temperatures. 
Structures that appear on the low energy side can not be seen with a free carrier model. 

This manifests a further example of the complicated interplay between the different levels 
in the QCL structure, which is not present in simple two-level systems. 

Another interesting remark is that, for larger dephasing, the counter-rotating terms, 
usually very small and ignored [„5], start to be relevant. The Coulomb shifts have opposite 
directions for counter rotating terms than for rotating wave contributions, and they should 
yield interesting effects in structures specifically constructed to study their relevance. 
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IV. SUMMARY 



In summary, our nonequilibrium many-body theory yields a surprising contrast between 
the importance of Coulomb effects at the absorption and gain region in quantum cascade 
structures. While these effects are small in the gain region for the structure considered, 
strong features appear at THz absorption. 

The many-body effects are influenced by both nonequilibrium occupation functions and 
Coulomb matrix elements, which depend on wavefunction overlap. Those two factors can 
be modified by means of an external bias which ultimately allows us to control the overall 
strength of Coulomb corrections in the structure. 
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APPENDIX A: INTERPOLATION FORMULA FOR COULOMB MATRIX EL- 
EMENTS 

The Coulomb matrix elements in Eq. |3] can be difficult to evaluate numerically, and 
we give here tractable expressions for the relevant terms. In general, the Coulomb matrix 
elements are given by 



Vl^^) = V?F( ti "_ X F), (AI) 



F[^ VX ^\ = I dz dz' <t>*Jz)Mz) ex P - A) <f>*x(z')Mz'), 



q j \ q 

q 

q ' e Aq- 

As pointed out in the text, Vq^ v matrix elements can be expressed exactly as a depo- 
larisation length in the small wavenumber limit, 



The large wavenumber limit reads 



V (»H">\ =V ~1 (A3 ) 



±- = 2 / dz\Mz)\ 2 \Mz)\ 2 - 

Inspection of the two extremes, suggests an interpolation formula for the depolarisation- 
like terms, 



Similar considerations yield an analogous interpolation formula for the subband shift and 
exchange matrix elements, 

v { V T) = V T 2 *-' < A5) 

V Q J 1 + k^q 

where 

L * = J^t£wF' (A6) 

L ex _ ' 



~" fdz\M*)\ 2 \M*)\ 2 ' 

The lengths, , Ll x , L h ^ v have been defined so that they reduce to the quantum well 
width if the wavefunctions used are the eigenstates of an infinite potential well. 
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TABLE I: Absorption/gain transition parameters. The first column is in Volts/period, (v,fi) de- 
notes the lower and higher subbands in a transition, n u — n^ is the difference in subband population 
in units 10 11 /cm 2 . Lf, x , L„ 9 U , and Lj^ are in nm and p UfJi in exnm. — E v is the energy 

difference in meV at k=0 of Eq. (4), i.e. including the static mean field corrections, but no further 
Coulomb effects. Y Vil = T u + T M is the dephasing in meV. 
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FIG. 1: Section of the quantum cascade laser structure and a few relevant levels within a period, 
including mean field effects due to the ionised dopants in the structure, and for a bias V=0.3 
V/period. Level 7 is an injector level. Levels 6 and 8 yield the main contribution for the midinfrared 
gain spectrum, while levels 7 and 9 are responsible for the main THz absorption feature. 
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FIG. 2: Influence of bias on absorption/gain spectra. From top to bottom the bias is V=0.3, 0.2 
and 0.12 V/period. In all panels, the dot-dashed lines are for free carriers, while the solid curves 
have full Hartree-Fock Coulomb corrections. 
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FIG. 3: Details of THz absorption including mean field effects in the dominant transition of each 
bias. From bottom to top, the bias is V= 0.12, 0.2 and 0.3 V/period, and the corresponding 
transitions are between levels (6,7), (7,8) and (7,9) as discussed in the text. In all panels, the solid, 
dotted, short-dashed, long-dashed and dot-dashed curves are respectively for Full Hartree-Fock, 
sub-band shift only, depolarisation only, exchange only, and free carriers. 
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FIG. 4: Details of the main transitions that contribute to the gain spectra including mean field 
effects for a bias V= 0.3 V/period. The lower and top panels are respectively for (6,8) and (6,7) 
transitions. In both panels, the solid, dotted, shot-dashed, long-dashed and dot-dashed curves are 
respectively for Full Hartree-Fock, sub-band shift only, depolarisation only, exchange only, and free 
carriers. 
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FIG. 5: Comparison of absorption spectra with (solid) and without (dot-dashed) Hartree-Fock 
Coulomb corrections. In all panels the bias is V=0.2 V/period. The carrier densities are given 
respectively by N=3.5, 4.6, 5.25, and 7 x 10 n /cm 2 . 
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FIG. 6: Gain spectra as a function of doping. All curves in the top panel have full Hartree-Fock 
corrections. The solid, dotted, long-dashed and dot-dashed curves correspond, respectively to 
N=3.5, 4.6, 6.4 and 7 x 10 11 /cm 2 . The inset depicts the spectral position of peak gain (in meV) 
as a function of doping in 10 11 /cm 2 . The open squares are for full Hartree-Fock while the circles 
are for free carriers. The symbols are connected by solid and dot-dashed lines respectively to guide 
the eye. In the central and bottom panels, the dopings are, respectively, N=3.5 and 7 x 10 11 /cm 2 . 
The solid and dot-dashed curves are respectively for full Hartree-Fock and free-carrier calculations. 
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FIG. 7: Comparison of absorption spectra with (solid) and without (dot-dashed) Hartree-Fock 
Coulomb corrections for a bias of 0.2 V/period . The top and lower panels are respectively for 
phonon bath temperatures T=77 and 180 K. 
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